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tj ' ABSTRACT 

o' 

Aims. We aim at reproducing the mass- and o--[afFe] relations in the stellar populations of early-type galaxies by means of a cosmo- 
logically motivated assembly history for the spheroids. 
CO ' Methods. We implement a detailed treatment for the chemical evolution of H, He, O and Fe in GallCS, a semi-analytical model for 

galaxy formation which successfully reproduces basic low- and high-redshift galaxy properties. The contribution of supemovae (both 
type la and II) as well as low- and intermediate-mass stars to chemical feedback are taken into account. The model predictions are 
compared to the most recent observational results. 

Results. We find that this chemically improved GallCS does not produce the observed mass- and cr-[a/Fe] relations. The slope is too 
^ , shallow and scatter too large, in particular in the low and intermediate mass range. The model shows significant improvement at the 

. highest masses and velocity dispersions, where the predicted [a/Fe] ratios are now marginally consistent with observed values. We 

' show that this result comes from the implementation of AGN (plus halo) quenching of the star formation in massive haloes. A thor- 

' ough exploration of the parameter space shows that the failure of reproducing the mass- and cr-[a/Fe] relations can partly be attributed 

to the way in which star formation and feedback are currently modelled. The merger process is responsible for a part of the scatter We 
suggest that the next generation of semi-analytical model should feature feedback (either stellar of from AGN) mechanisms linked to 
single galaxies and not only to the halo, especially in the low and intermediate mass range. 
^ ' Conclusions. The integral star formation history of a single galaxy determines its final stellar [ff/Fe] as it might be expected from the 

CO , results of closed box chemical evolution models. However, the presence of dry-mergers and metal recycling in the hot gas phase helps 

■ in keeping the a element abundance in the stars at a super-solar level in a hierarchical galaxy formation scenario. 

!>■ ■ 
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1. Introduction predicted star formation histories of massive elliptical galaxies 

were to extended to produce a/Fe consistent with observations. 

The Cold Dark Matter (CDM) scenario (Peebles, 1982) success- This conclusion was later reinforced by Nagashima et al. (2005), 

, ^ fully explains the growth of the large scale structure of the uni- who included a self-consistent treatment of chemical enrich- 

k> verse (Springel, Frank & White 2006). Since the original appli- ment in semi-analytic galaxy formation models, and by Pipino 

j_j cation of this scenario to the galaxy formation process problem & Matteucci (2006, PM06). 

C3 (Kauffmann & White, 1993, Cole et al 1994), however, several More recendy, a plethora of new models have been presented 

modifications to the pure hierarchical assembly of the building to address the general phenomenon of downsizing (e.g. Croton 

blocks had to be introduced in order to deal with the complex- et al 2006, De Lucia et al. 2006, Bower et al 2006, Cattaneo et 

ity of baryonic physics. Among the main open issues, we men- al., 2005, Sommerville et al. 2008, Kaviraj et al. 2005, Fontanot 

tion the anti-hierarchical behaviour of the AGNs (e.g. Hasinger et al., 2007). The new key ingredient in these models is feedback 

et al. 2005), the evolution of luminosity function with red-shift from super-massive black holes, which is used to suppress resid- 

(e.g. Bundy et al. 2005) as well as the increase of mean stel- ual star formation at late times in the evolution of massive galax- 

lar [a/Fe] with galaxy mass (or cr) in elliptical galaxies (e.g., ies (Granato et al. 2004). Such a suggested scenario seems to be 

Worthey et al 1992; Trager et al 2000; Thomas et al 2005, Nelan supported by observations (e.g. Nesvadba et al. 2006). In prac- 

et al. 2005).This relationship together with the old infeiTed ages tice, the mass assembly still occurs at late times in these models 

implies that more massive ellipticals formed earlier and faster ^ whereas most of the stars have been formed at high red-shift 

with respect to smaller objects (Matteucci 1994; Thomas et al. small sub-units (but see Cimatti et al. 2006). The preferred 

2005). In the following we will refer to the above observational mechanism for the assembly of massive spheroids is a sequence 

constt-aints with the term downsizing (Cowie et al. 1996). of dissipation-less (dry) mergers. This leads generally to better 

Thomas (1999) was the first to study the chemical enrich- agreement with the observed downsizing pattern, 

ment of a and Fe-peak elements in the framework of hierarchi- However, the [a/Fel-mass relation has not been studied in 

cal models of galaxy formation. In a very simplistic approach, these new generation models (see Pipino & Matteucci 2008, 

Thomas (1999)ran chemical evolution simulations over the star PM08). The aim of this paper is to fill this gap. To this end we 

formation histories predicted by Kauffmann (1996) neglecting implement a fully self-consistent treatment of the chemical evo- 

the complex merger history of the galaxies. It turned out that the lution, which includes a robust estimate of the type la supernova 



A. Pipino et al.: GallCS II: [a/Fe] ratios in eUipticals 



2 

rate and of the Fe production into GallCS (Hatton et al., 2003; 
Paper I hereafter), a state-of-the-art semi-analytical model for 
galaxy formation and evolution based on a CDM-driven growth 
of the structures. The main goal is to check model predictions 
against the latest observational results for the a/Fe-mass rela- 
tion and the mass-metallicity relation (MMR, e.g. Carollo et al. 
1993), that so far could be simultaneously accounted for only by 
revised monolithic models (Pipino & Matteucci, 2004, PM04). 
The new results will be then interpreted in the light of our pre- 
vious work with the chemical evolution models (Thomas 1999, 
PM04, PM06, PM08). 

The structure of the paper is as follows: the main improve- 
ments with respect to Paper I are spread out in Sec.2; in Sec. 3 the 
chemical evolution scheme is tested against the Milky Way and 
the local SNIa rate. In Sec. 4, 5 and 6 the results are presented 
and discussed, respectively. 

2. The model 

2.1. The GallCS galaxy formation model 

GallCS is a model of hierarchical galaxy formation which 
combines high resolution cosmological N-body simulations to 
describe the dark matter content of the Universe with semi- 
analytic prescriptions to follow the physics of the baryonic mat- 
ter. GallCS has been thoroughly presented in Paper I, which we 
refer the reader to for a detailed discussion of the model assump- 
tions and properties. It has already been used for the study of 
the colour-magnitude relation and the progenitor bias of ellipti- 
cal galaxies (Kaviraj et al. 2005), the reproduction of the Galex 
NUV-optical colours (Kaviraj et al. 2007) and the black-hole 
mass - cr relation (Cattaneo et al 2005). It has also been used to 
explore the consequences the halo-quenching mechanism (Keres 
et al. 2005) by Cattaneo et al. (2008 and references therein). The 
above mentioned papers represent a comprehensive set of bench- 
mark tests that we will not repeat here, but simply point out that 
results are preserved to a large extent in our present implemen- 
tation. 

We briefly recall the specifications of cosmological N-body 
simulation used to construct the halo merger trees. This simu- 
lation is a realization of flat cold dark matter universe with a 
cosmological constant of Qa = 0.667. The simulated volume is 
a cube of side Lbox = IOO/ZkJqMpc, with him = //o/lOOkm s"' - 
0.667, which contains 256"* particles of mass 8.3 x lO'^M© each, 
the cold dark matter power spectrum was normalised in agree- 
ment with the present day abundance of rich clusters (erg = 
0.88). One should bear in mind that the dark matter simulation 
cannot resolve haloes less massive than 1.6 x IO^'M©, which 
implies that a galaxy less massive than 2 x lO^^M© is formally 
below the resolution limit. The spatial resolution, instead, is such 
that we cannot resolve scales below 30 kpc. 

As hot gas cools and falls to the centre of its dark matter halo, 
it settles in a rotationally supported disc. According to Paper I, if 
the specific angular momentum of the accreted gas is conserved 
and starts off with the specific angular momentum of the dark 
matter halo, we assume it forms an exponential disc with scale 
length rd given by: 

A 

rd = —pRioQ- (1) 

Galaxies remain pure discs if their disc is globally stable (i.e. 
Vc < 0.7 X V,oi where V,,,, is the circular velocity of the disk- 
bulge-halo system ; see e.g. van den Bosh et al 1998), and they 



do not undergo a merger with another galaxy. In the case where 
the latter of these two events occurs, we employ a recipe to dis- 
tribute the stars and gas in the galaxy between three components 
in the resulting, post-merger galaxy, that is the disc, the bulge, 
and a star-burst (see Paper I). In the case of a disc instability, we 
simply transfer the mass of gas and stars necessary to make the 
disc stable to the burst component, and compute the properties of 
the bulge/burst in a similar fashion as that described in Paper I. 
Bulges are assumed to have a density profile given by Hernquist 
(1990). The bulges are assumed to be pressure supported with a 
characteristic velocity dispersion cr, computed at their half-mass 
radius. 

Galaxy morphology in the model is determined by the ratio 
of the B-band luminosities of the disc and bulge components. A 
morphology index is defined as 

/ = exp(^) (2) 

such that a pure disc has 1=1 and a pure bulge has 7 = 0. 
Following Baugh et al (1996), ellipticals have I < 0.219, SOs 
have 0.219 < I < 0.507 and spirals have I > 0.507. This sim- 
ple prescription is clearly incapable of capturing the complex 
spectrum of real galaxy morphologies. Therefore, in what fol- 
lows, 'spirals' refer to all systems which do not have a dominant 
spheroidal (bulge) component. ObservationaUy, this includes not 
only systems with distinctive spiral morphologies, but also pecu- 
liar or irregular systems. 

2.2. The chemical evolution 

The main novelty of the present versions of GallCS is the im- 
plementation of a self-consistent treatment of the chemical evo- 
lution with finite stellar lifetimes and both type la and type II 
supernovae ejecta. In practice, we follow the chemical evolu- 
tion of only four elements, namely H, He, O and Fe. This set 
of elements is good enough to characterise our simulated ellipti- 
cal galaxy from the chemical evolution point of view as well as 
small enough in order to minimise computational resources. In 
fact, as shown by the time-delay model (Matteucci & Greggio, 
1986), the [a/Fe] ratio is a powerful estimator of the duration 
of the SF. Moreover, both the predicted [Fe/H]-mass and [Z/H]- 
mass relationships in the stars can be tested against the observed 
Colour-Magnitude Relations (hereafter CMRs; e.g. Bower et al. 
1992, Kaviraj et al 2005) and MMR. In order to clarify this 
point, we recall that the O is the major contributor to the total 
metallicity, therefore its abundance is a good tracer of the metal 
abundance Z. Moreover, in this paper we focus on the theoretical 
evolution of the a elements, and the O is by far the most impor- 
tant. On the other hand, the Fe abundance is probably the most 
commonly used probe of the metal content in stars, therefore it 
enables a quick comparison between our model predictions and 
the existing literature. 

In the foUowing [a/Fe] ratio will always refer to the 
luminosity-weighted average over the stellar populations that 
make a galaxy, unless stated otherwise. This value guarantees 
a robust comparison with its observational counterpart, namely 
to the "SSP-equivalent" value inferred from the integrated spec- 
tra of elUptical galaxies. We refer to Pipino et al. (2006) for de- 
tails and caveats on the use of SSP-equivalent" abundances and 
abundance ratios as proxies for the mean properties of a com- 
posite stellar populations like an elliptical galaxy. 
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The star formation rate in the disc is 



iA(f) = 



(3) 



yn 



Here Mcoid is the mass of the gas in the disc (all the gas in the disc 
is cold and all the gas in the halo is hot) and fdyn is the dynamical 
time (the time to complete a half rotation at the disc half mass 
radius). The parameter yS,, which determines the efficiency of 
star formation has a fiducial value of /J, = 50 (Guiderdoni et al 
1998). 

A Salpeter (1955) initial mass function (IMF) constant in 
time in the range 0.1 - AQMq is assumed, since PM04 showed 
that the majority of the photochemical properties of an ellip- 
tical galaxy can be reproduced with this choice for the IMF. 
We adopted the yields from Iwamoto et al. (1999, and refer- 
ences therein) for both SNIa and SNII. The SNIa rate for a SSP 
formed at a given radius is calculated assuming the single degen- 
erate scenario and the Matteucci & Recchi (2001) Delay Time 
Distribution (DTD). The convolution of this DTD with if/ (see 
Greggio 2005) gives the total SNIa rate, according to the follow- 
ing equation: 



Rla(t) 



c 



■min(t,T:^) 



A{t - T)i//{t - T)DTD{T)dT 



(4) 



where A{t-r) is the fraction of binary systems which give rise to 
Type la SNe. Here we will assume it constant (see Matteucci et 
al. 2006 for a more detailed discussion). The time r is the delay 
time defined in the range (t,, t J so that: 



r 



DTD(T)dT = 1 



(5) 



where t, is the minimum delay time for the occurrence of Type 
la SNe, in other words the time at which the first SNe la start 
occurring. We assume, for this new formulation of the SNIa rate 
that T, is the lifetime of a 8M0 star, while for Tj, which is the 
maximum delay time, we assume the lifetime of a O.8M0 star 
Finally, ka is the number of stars per unit mass in a stellar gen- 
eration and contains the IMF. The detailed treatment of SNIa is 
a substantial improvement with respect to Paper I. 

Stars are evolved between time-steps using a sub-stepping of 
at least 1 Myr During each sub-step, stars release mass and en- 
ergy into the interstellar medium. In GallCS, the enriched mate- 
rial released in the late stages of stellar evolution is mixed to the 
cold phase, while the energy released from supernovae is used 
to re-heat the cold gas and return it to the hot phase in halo. The 
re-heated gas can also be ejected from the halo if the potential 
is shallow enough (see also Paper I). The rate of mass loss in 
the supernova-driven wind that flows out of the disc is directly 
proportional to the supernova rate. 

The original formula for the chemical processing of the total 
metal content (see Paper I) has been extended to the elemental 
species we deal with, so that the ejecta in the gas mass from the 
stellar population are: 

£/(f) - I t//(t-t„)([m-w(m)]Zi^coidit-tm)+mYi(m))<p(m)dm(6) 

Jm{l) 

where m{t) is the mass of a star having lifetime w{m) is the 
mass of the remnant left after the star has died, and (f){m) is 
the IMF. The first term on the right hand side represents the 
re-introduction of the metals that were originally in the stars 
when they formed, and y,(m) is the fraction of the initial stel- 
lar mass transformed via stellar nucleosynthesis into the element 



i. Throughout this work, we assume chemical homogeneity (in- 
stantaneous mixing), such that outflows caused by feedback pro- 
cesses are assumed to have the same metallicity as the inter- 
stellar medium, though in reality the material in the outflow is 
likely to be metal-enhanced (see Sec 6.3). 



2.3. Galaxy evolution and properties 

The fundamental assumption is that all galaxies are born as discs 
at the centre of a dark matter halo. The transformation of disc 
stars into bulge stars and of disc gas into star-bursting gas is due 
to bar instabilities and mergers. Gas is never added to bulges 
directly and the only gas in bulges is that coming from stellar 
mass loss. The star-bursting gas forms a young stellar population 
that becomes part of the bulge stellar population when the stars 
have reached an age of 100 Myr. We do not readjust the bulge 
radius when this happens. The disc has an exponential profile, 
while the bulge and the star-burst are described by a Hernquist 
(1990) density distribution. The star-burst scale is rburst = ''''"bulge 
with K = O.l. 

The star formation law (Eq. [3]) has the same form and uses 
the same efficiency parameter /3, for all three components when 
we redefine Mcoid as the mass of the gas in the component and 
fdyn as the dynamical time of the component. For the components 
described by a Hernquist profile, the dynamical time is fdyn = 
ro s/cr, where ro.5 is the half mass radius and cr is the velocity 
dispersion at the half mass radius. 

The fraction of the disc mass transferred to the spheroidal 
component (the bulge and the star-burst) depends on the mass 
ratio of the merging galaxies. The separation between a minor 
and a major merger is for a mass ratio of 1:3. 

GallCS cannot spatially resolve galaxies, therefore we can 
only predict chemical properties averaged over the galactic ra- 
dius. 



2.4. Energetics 

The SNII feedback is given by: 



m = 2i//(t) T 



(7) 



where e is the efficiency of the supernova-triggered wind which 
is proportional to v^,^. and depends both on the porosity of the 
ISM (see Silk 2001 for details) and the mass-loading factor This 
latter accounts for entrainment of interstellar gas by the wind and 
can be considered as a free parameter whose value is around 10. 
Note that in the previous equation, t/sn is the number of super- 
novae per unit star-forming mass, which is a prediction of the 
Initial Mass Function (IMF) chosen, and Esn is the energy of a 
supernova, assumed to be 10^'erg. 

At variance with chemical evolution models as PM04, where 
the total la-HlI SNe feedback is sufficient to halt the SF, Paper I 
relies onto the observed correlation between AGN and velocity 
dispersion (Ferrarese & Merrit, 2000), and simply prevents gas 
from cooling in a halo which as a mass above the critical value of 
~ 10" M0 to quench cold gas accretion (Granato et al. 2004). A 
further halo-quenching mechanism has been implemented into 
GallCS by Cattaneo et al. (2008), who showed how this further 
refinement leads to a better agreement between our model pre- 
dictions and SDSS observations of the luminosity function and 
the colour bimodality (Baldry et al. 2004). 
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Fig. 1. Dotted lines: metallicity distribution of stars as functions 
of [Fe/H] and [O/Fe] (3cr contour) obtained by stacking the MW- 
like spirals predicted by GallCS. Solid line: [O/Fe] as a function 
of [Fe/H] as predicted by GallCS for a median MW-like spiral. 
Dots: observational data in the range of interest as compiled and 
homogenised by Frangois et al. (2004). 



3. Calibration of the model 

3. 1 . Comparison to the Milky Way 

To provide a consistency check for our model, we adopt a pro- 
cedure typical of chemical evolution studies. We first compare 
MW-like galaxies in our simulations with the known properties 
of our galaxy. Since many properties of such galaxies have al- 
ready been tested in the calibration of Paper I, here we use the 
same selection criteria of MW-like galaxies (namely mgas/ffJbar - 
0.10+0.05, = -23.7mag±0.3mag,ye = 220 km/ s +20 km/ s 
and requiring that the galaxy have spiral morphology) and we 
show only chemical evolution predictions. We found that ~3% 
of the spiral population is made by MW-like objects in agree- 
ment with Paper I statistics. 

In Fig. [U we plot [O/Fe] ratio as a function of [Fe/H] in 
the stars of the MW-like spirals predicted by GallCS. The dot- 
ted lines give the distribution of stars formed out of gas with 
a given chemical pattern (i.e. a given {Fe/H~\ and {0/Fe\) as 
a 3cr contour in the [O / F e\-[F e / H] plane, whereas the thick 
solid line give the median trend. The agreement with the over- 
all trend observed in our own Galaxy makes us confident that 
the model is correctly calibrated. Unfortunately, due to the 
metallicity resolution (five bins in total metallicity Z, namely 
0.001,0.004,0.008,0.02 and 0.04) of the code and to the fact that 
galaxies are identified only once their host DM haloes are quite 
massive (M,,,> > 1.6 x IO'^Mq), we cannot explore the region 
at [Fe/H] < -1.5. We plan to use higher mass resolution simu- 
lations in future work, but for sake of comparison with Paper I 
we restrict our analysis to the same simulation that was used in 
Hatton et al. (2003). 




11.0 11.5 

Log stellar mass 



Fig. 2. Present-day SNIa rate in SNuM units as a function of 
the galactic mass for our model ellipticals (black crosses). The 
hatched are brackets the Icr-scatter region around the observa- 
tional fit (0.066 SNuM, Mannucci et al., 2008). 



3.2. Present-day SNIa rate 

We then verify that the model galaxies predict a present-day 
morphology dependent SNIa rate in agreement with observa- 
tions. The MW-like spirals presented above exhibit a SNIa rate 
of 0.09 SNuM (i.e. specific SN explosion rate in units of IO'^Mq 
of stars per century), which is in fair agreement with the obser- 
vational estimates (0.06!5J{J['^ for SOa/b and 0.14+°°45 ^j. ^^^.^^^ 
respectively, see Mannucci et al. 2008) given the fact the GallCS 
does not allow a finer morphological classification. As shown 
in Fig. |2] the vast majority of our simulated ellipticals exhibit a 
present-day SNIa rate within Icr from the observational mean 
value given by Mannucci et al. (2008). This test basically guar- 
antees that, given the star formation history of the model galax- 
ies, we have calibrated the uncertainties in the progenitor nature 
and delay time distribution of SNIa which are incorporated in the 
parameter A (see Eq.|2]i. In particular, in order to reproduce the 
present-day observed SNIa rate we assume A - 0.0025 which 
is the value typically adopted in chemical evolution models of 
the Milky Way (see Matteucci et al. 2006). As a result the Fe 
production rate from SNIa is also calibrated. 



4. The cr- and mass-[Qf/Fe] relations 

4.1. The standard GallCS model 

From this section onwards we deal with the mean novelty of the 
present work, namely the study of the predicted a/Fe-mass rela- 
tion and its comparison to the observations. To be consistent with 
observed values we present luminosity-weighted values which 
take into account the disc component (if any). We stress, how- 
ever, that the mass-weighted quantities do not differ much from 
the luminosity-weighted ones especially at the high mass end of 
the sample, where SF has been suppressed at high redshift. At 
variance with our previous work (PM04, PM06, PM08) on the 
ff-enhancement in the very central part of the galaxies, here we 
present our predictions on the [a/Fe] ratio in the whole galaxy 
(we recall that GallCS's spatial resolution is larger than typi- 
cal galactic effective radii) and consistently compare them to the 
recent observational estimates by Thomas et al. (2008) which 
pertain to the entire galaxies. 

We also recall that observations suggest that the observed ra- 
dial gradient slope in the [(r/Fe] has, on average, a null value (e.g. 
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Fig. 3. The a/Fe-mass and -cr relations as predicted by GallCS 
for the whole sample of ellipticals (black points). The thick solid 
lines encompass the Icr-region (hatched area) around the mean 
trend reported by Thomas et al. (2008). 
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Fig. 4. The a/Fe-mass relation as predicted by GallCS for the 
whole sample of ellipticals (black points). A subsample of ellip- 
ticals older than 10 Gyr is emphasised with a lighter colour. The 
thick solid lines encompass the Icr-region (hatched area) around 
the mean trend reported by Thomas et al. (2008). 



Mehlert et al. 2003). As shown by Pipino, D'Ercole & Matteucci 
(2008a), in fact, even though the observed gradients (CaroUo et 
al. 1993, Davies et al. 1993) suggest that most ellipticals form 
outside-in, the expected strong and positive [a/Fe] gradient can 
be affected by the metal rich gaseous flows inside the galaxy 
acting together with the SFR. The net result is a gradient in the 
[a/Fe] ratio nearly flat. Hence, we can safely neglect the pres- 
ence of gradients in our study. Instead, they might affect the 
MMR (see discussion in Sec.|5]i 

The results for our fiducial GallCS version are presented in 
Fig. [3] In agreement with Nagashima et al. (2005), the a/Fe ra- 
tios do not show any correlation with mass (Fig. [3] top panel), in 
strong contrast with the clear positive correlation derived obser- 
vationally. The simulations seem to produce decreasing [a/Fe] 
ratios with increasing galaxy mass, with a slight upturn at the 
high-mass end. As a result the scatter is large. A linear fit of the 
simulation results in the [a/Fe]-mass plane would give a flat rela- 
tionship. It is interesting to note that this result does not depend 
on the environment in the sense that, if we restrict the regres- 
sion analysis to a sub-sample of galaxies living in haloes whose 
mass is comparable to rich cluster of galaxies, we do not notice 
substantial changes in the predicted relationships. 

We notice that the most massive galaxies attain a typical level 
of a-enhancement that is only Icr off the value suggested by the 
observations. This is an improvements with respect to previous 
results and mainly caused by the implementation of AGN feed- 
back (see Discussion). There are two possible formation paths 
for these objects: i) either these galaxies assemble through dry 
(gas-poor) mergers or ii) assemble most of their mass over very 
short time-scales (less than 0.5 Gyr). Indeed, this ensures that 
the pollution from SNIa is kept at a low level and, hence, that 
they maintain an over-solar [a/Fe] ratio in their stars. However, 
as showed by Pipino & Matteucci (2008), low-mass and highly 
a-enhanced galaxies are needed if one wants to create the most 
massive spheroids with a suitable a enhancement by means of 
dry-mergers. GallCS then predicts that a small number of them, 
with masses ~ 0.5 - 1 x IO^Mq should survive down to red- 
shift zero. Unfortunately such galaxies are not observed. In a 
sense one could turn the argument around and say that a robust 
prediction of semi-analytic models of hierarchical galaxy forma- 
tion is the presence of low-mass, highly [a/Fe] -enhanced galax- 
ies at high redshift because it is the only way in these models 
to build local massive ellipticals with the observed [a/Fe] ratios. 
We know from observations that such objects do not seem to ex- 
ist at moderately high redshifts around z ~ 0.4 (Ziegler et al. 
2005) 

If we consider the subsample of galaxies whose luminosity- 
weighted ages are larger than 10 Gyr (lighter points in fig.llli we 
notice that the galaxies populating the region below the observed 
area in the [a/Fe] -mass plane disappear 

We attribute this to the fact that these galaxies live in the 
centre of massive haloes where both the original AGN feed- 
back implemented in GallCS and the halo-quenching mecha- 
nism (Cattaneo et al. 2008) halted the cooling in the gas a long 
time ago. This finding confirms the interpretation of the [a/Fe] 
ratios below the observed range as being related to a too long du- 
ration of the star formation, as we will discuss later in the paper. 
However, there has been some improvement in the agreement 
with observations with respect to previous models. 

Similar results are obtained when plotting the [a/Fe] as 
a function of the stellar velocity dispersion cr (fig. [3] bottom 
panel). Comparing the two panels in fig. [3] we notice that the 
scatter is somehow reduced and that the galaxies follow a trend 
which is closer to the observational results. We can understand 
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Fig.5. Integral star formation history for a ~ lO^M© galaxy 
with strong a-enhancement (solid line). The SFH inherited from 
the single building blocks is also shown by dotted (progenitors 
merging very early on), dashed (progenitors merging at z~4) and 
long dashed (progenitors merging at z~2-3) lines, respectively. 



this as we expect the velocity dispersion to be more correlated 
with the properties of the DM host haloes, whereas the baryonic 
mass is more sensitive to our modelling of feedback processes. 
However, one should bear in mind that GallCS assumes virial- 
ization and a fixed density profile (see Paper I) to calculate cr. 
Whilst this assumptions are reasonable at z=0, they very likely 
should be revised at high redshifts. 

Moreover, we should note that the scatter is still much larger 
than the observed one and that, for a given velocity dispersion, 
the model galaxies tend to have on average a lower \ajFe\ ra- 
tio than the observed ones. The latter problem can be handled 
in several ways (stellar yields, IMF, feedback, SF efficiency, see 
Sec.|6]for details), whereas the former is intrinsic to the model, 
being linked to the stochastic nature of the merger process where 
errors in the estimate of individual SF histories add up as galax- 
ies merge together. 

In order to better understand the origin of this scatter, we 
now focus on few selected galaxies with the same mass, but very 
different [O/Fe] and we study their star formation histories. We 
perform this exercise for two masses: i) typical ellipticals of ~ 
IO^Mq; ii) massive ellipticals of ~ 6 x IO'^Mq. 

4.2. Typical ellipticals 

We start the analysis in the ~ 10" M© mass range, where most of 
the predicted galaxies scatter outside of region of the observed 
values. A useful test bench for understanding the behaviour of 
such galaxies is provided by running chemical evolution mod- 
els with the same stellar yields, IMF and the same integral star 
formation rate as the selected semi-analytic galaxies. 

Indeed when we force the SFH pictured in Fig.|5](solid line) 
to happen in a standard chemical evolution model (PM04), the 
predicted [O/Fe] is lower by 0.2 dex than that obtained with 



Fig. 6. Integral star formation history for a ~ 10"Mq galaxy 
with nearly solar [a/Fe] 
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Fig. 7. The a/Fe-mass relation as predicted by GallCS for the 
whole sample of ellipticals (black points). A subsample of ellip- 
ticals which experienced at least two mergers is presented with 
lighter colours. The thick solid lines encompass the Icr-region 
(hatched area) around the mean trend reported by Thomas et al. 
(2008). 

GallCS. This brings the galaxy from a value of [a/FeJ^O.Sl 
down to a value of [a/Fe]=0.12 which is within the range of 
the observed values. 

On the other hand, for the SFH presented in Fig.|6]we have an 
ff-depletion with GallCS, which does not occur in PM04. Once 
again this change is enough to bring the galaxy back within the 
range of observed values. 

The difference here is that the galaxy whose SFH is por- 
trayed in Fig. |5] has 9 progenitors with 4 of them merging very 
early on (i.e. at z < 4.7), and which is passively evolving from 
red-shift 2. Looking at Fig. |5] (dashed and dotted lines) we see 
that all the progenitors of this galaxy actually have individual SF 
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time-scales which are shorter than the one that would be inferred 
from the mass-weighted SFH of the galaxy itself. This explains 
why the [or/Fe] ratio calculated by GallCS is higher than the one 
derived by a pure chemical evolution model of a single object 
with the same mass-weighted SFH. This is a systematic trend of 
hierarchical galaxy formation when compared to pure chemical 
evolution models as is shown also in Figs. |8]and|9]for massive 
ellipticals. On the other hand, in the case of fig. |6] we have only 
one progenitor which explains why the [a/Fe] ratios of GallCS 
and the pure chemical evolution model are in good agreement. 
Looking at statistics with the help of Fig.|7]the latter case (black 
points) represents 42% of the total number of elliptical galaxies 
and it is biased towards lower mass as we would expect since 
massive ellipticals are built by multiple mergers in the hierarchi- 
cal galaxy formation scenario (lighter points in Fig. |7] indicate 
galaxies which experienced at least two mergers.). One might be 
worried that these results depend on the mass resolution of the N- 
body simulation, however as stated earlier in the paper our mass 
resolution is such that galaxies more massive than 2 x IO^'Mq 
are resolved and this mass is about a factor of 10 lower than the 
mass of ellipticals considered in our analysis. This means that 
increasing the resolution would only change the number of mi- 
nor mergers (except, of course, in the very early stages of the 
formation of these galaxies). 

4.3. Massive ellipticals 

If we analyse the most massive galaxies, we confirm the findings 
of the previous section (see figs. |8}|9]l. 

In general, despite the quite high a-enhancement, the inte- 
gral star formation history appears to be broader than the one 
obtained by means of a monolithic collapse model which would 
predict the same [a/Fe] ratio. To better explain this point, let 
us think to the ideal case in which we have two progenitors 
with the same masses, same star formation histories (such that 
their final [a/Fe] is appropriate for their mass), but say also that 



their peak in the star formation rate is shifted by about 1 Gyr 
(one is younger than the other and this difference in age cannot 
be detected with the standard line-strength indices technique if 
this objects are more than 10 Gyr old). Let us also assume that 
these two galaxies coalesce via a dry merger later in their evolu- 
tion. The final object still has the same [a/Fe] of the progenitors, 
its mass in only doubled (therefore it still matches the observa- 
tions, given the spread in the a/Fe-mass relation, see PM08), 
but the integral star formation history will look like the one in 
Fig. |8] thus broader than the 0.5-0.7 Gyr expected by PM04 
(their model II). 

We also note that in GallCS, the galaxies do not evolve as 
closed box. They instead exchange metals with the surrounding 
hot halo as well as stars can be created in discs and the moved to 
bulges because of instabilities. Such processes, render the inter- 
pretation of the final [a/Fe] ratio on the basis of the SFH alone 
much more complicated. 

Finally, a stacked specific (i.e. per unit stellar mass) SFH for 
galaxies in different mass bins has already been presented by 
Cattaneo et al. (2008) and we do not repeat the analysis here. 
With the help of the SFH presented in this section, however, we 
explained why the average duration of the SF is a factor of 3- 
5 longer (and consequently the peak value is at a factor of 3-5 
lower) than what is required from pure chemical evolution stud- 
ies on line-strength indices analysis to reproduce the [a/Fe] in 
massive ellipticals. 

5. The mass-metallicity relation and ages of the 
galaxies in the standard model 

In fig. [To] we show the predicted MMR relation against the data 
by Thomas et al. (2008). We note that the predicted flat be- 
haviour is unaff'ected by the cuts in either mass or age as done in 
the SecH 

A failure in reproducing the MMR is expected on the basis 
of the preliminary analysis by PM08 who showed that the diag- 
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Fig. 10. The MMR as predicted by GallCS for the whole sam- 
ple of ellipticals (black points). The thick solid lines encompass 
the Icr-region (hatched area) around the mean trend reported 
by Thomas et al. (2008). Note that elliptical galaxies do ex- 
hibit quite strong [Z/H] gradients within one effective radius (e.g 
Carollo et al., 1993; Davies et al 1993) and Thomas et al. (2008) 
galaxies were observed with a fixed fiber size. Therefore it is dif- 
ficult to make a meaningful comparison between our predictions 
and observations as in the case of the [a/Fe]-mass relation (see 
text). 



nostic power of the MMR and a/Fe-mass relation relies in the 
fact that the mechanisms required to satisfy the former tend to 
worsen the agreement with the latter, and vice-versa. Only when 
most of the star formation process and the galactic assembly oc- 
cur at roughly the same time and the same place both relations 
can be fulfilled. In the past, in fact, models were mainly aimed in 
reproducing the MMR, rather than the downsizing trend, there- 
fore they failed in the a/Fe-mass relation. 

However, elliptical galaxies do exhibit quite strong (i.e. -0.3 
dex per decade in radius) [Z/H] and [Fe/H] gradients within one 
effective radius (e.g. Carollo et al., 1993; Davies et al 1993). 
Therefore it is difficult to make a thorough comparison between 
our predictions and observations. In practice, given the fixed 
aperture set by the SDSS fiber size, it is likely that smaller galax- 
ies contributed with most of their light, whereas only the cen- 
tral regions (more metal rich) are observed in bigger galaxies, 
thus biasing the observed slope of the MMR towards steeper 
values than the reality. However, a MMR does exist for the cen- 
tral regions of elliptical galaxies (e.g. Thomas et al. 2005) with 
quite a similar slope to the one holding for the entire galaxies. 
Therefore, the disagreement with the observed MMR at high 
masses cannot entirely explained by the fact that the GallCS spa- 
tial resolution is not enough to take the aperture effects into ac- 
count. We recall that, as shown by Pipino et al. (2008a, see also 
Sec. 12.31 ), this is not a problem for [a/Fe] ratios in the stars, be- 
cause of the interplay between internal metal flows and the star 
formation efficiency, which keep the [o'/Fe] gradient flat. 

We do not show predictions on the age-mass and colour mag- 
nitude relation. We refer the reader to Cattaneo et al. (2008) 
and Kaviraj et al. (2005) who show a remarkable agreement of 
the predictions made be means of GallCS and the latest obser- 
vational results. We only note that the predicted scatter in the 
predicted MMR and (SSP-equivalent) age-mass relationships is 
comparable to the intrinsic scatter derived by Thomas et al. 
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Fig. 11. a-enhancement versus number total number of mergers 
(top panel) and number of major mergers (bottom panel) for our 
model galaxies. 

(2005, 2008). On the other hand, the scatter in the a/Fe-mass 
relation is about twice as big. 

6. Discussion 

In the previous sections we showed the sensitivity of the [a/Fe] 
to the integral star formation history of the galaxy and estimated 
the offset in the predicted value with respect to pure monolithic 
formation. Moreover, we studied the scatter in the predicted 
[a/Fe] -mass relation in different mass ranges in order to under- 
stand the reasons for the persisting disagreement with observa- 
tions. In this section we further discuss this issue and try to find 
viable solutions to address it. 

Interestingly, according to fig.[TTl on average the final [a/Fe] 
ratio seems to be independent from the the total number of merg- 
ers (top panel) as well as from the number of major mergers 
(bottom panel). Moreover such a relationship also holds if we 
consider the number of either dry or wet mergers, and it is inde- 
pendent from the merger classificatiorQ. The model predicts a lot 

' Since gas is always present in the model galaxies, we tested three 
wet merger definitions - namely when the total mass in gas divided the 
total baryon mass of a newly formed galaxy exceed 0.01, 0.05 and 0.1 
- and we did not find substantial difi'erence among them. Note that dry- 
mergers are then defined by subtracting the number of wet mergers to 
the total for a given galaxy. 
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Fig. 12. ff-enhancement versus the epoch of the formation of the 
galaxies as given by the time at which the galaxy is spotted by 
GallCS, rather than from the SSP luminosity-weighted age. 

of 'monolithically' (i.e. mergers) behaving galaxies which ex- 
hibit a large scatter in the final [c/Fe] ratio. A closer inspection 
tell us that this fact has several reasons: 

i: Galaxies which evolve from disc to bulge morphology 
through instabilities at late times. 

ii: Galaxies in environments where the AGN-quenching thresh- 
old occurs too late. These galaxies (along with those of case 
i) populate the [a/Fe]-mass plane at values lower than the 
observed ones. In this case a better treatment of the feedback 
at galactic scales, possibly including (SNIa-driven) winds, 
might represent a solution. 

iii: Galaxies in environments where the AGN-quenching thresh- 
old occurs too early (these galaxies populate the [a/Fe] -mass 
plane at values higher than the observed ones). This issue 
seems to be common to other models of galaxy formation 
(Kimm et al., 2008). Since these galaxies are satellite of mas- 
sive haloes, one interesting (but probably not unique) possi- 
bility is to have them accrete some fresh gas over a long time- 
scale to keep SF going and, therefore, decrease the [a/Fe] 
ratio. 

There is a trend such that galaxies with more than ~ 5 merg- 
ers in their formation histories always have relatively high a/Fe 
ratios. As the most massive galaxies are those which underwent 
the highest number of both total and wet mergers in their life- 
time, these objects display relatively high a/Fe ratios and old 
ages. If we select again the galaxies with SSP-weighted age 
larger than 10 Gyr, we find that all the galaxies with a number 
of mergers larger than 7 belong to this category. (Not) surpris- 
ingly they are also the most massive (stellar masses larger than 
5 X 10"Mo). To be more precise, we now make use of the epoch 
of the actual birth (defined as the time at which the galaxy is first 
spotted by GallCS) of the galaxy, rather than the SSP-weighted 
age. Fig. [12] makes clear that the higher the redshift of the final 
assembly (estimated as the time at which the last merger occurs), 
the higher the final [a/Fe] ratio, as expected from the linear re- 
gression analysis on the predicted a/Fe-mass relation presented 
earlier. The majority of galaxies with a number of mergers in 
their lifetime lower than 4 (2) have formation epoch larger than 
1 (2) Gyr, and have mostly stellar masses smaller than 2 (1) 
xIO^'Mq. a confirmation of these findings comes from the pure 
chemical evolution study on gaseous mergers by PM06. They 



found that if a star-burst triggered by a significant accretion of 
pristine gas (comparable with the mass of stars already formed 
- roughly similar to a major wet merger in a galaxy formation 
picture) occuiTed at a significantly high redshift and just after 
the main burst of SF, the present day photo-chemical proper- 
ties of the final elliptical galaxy match the observed ones. Star- 
burst triggered at very low-redshift, instead, implied very low 
predicted final [a/Fe] ratios, even if the gas mass involved in the 
second burst was much lower than the one converted into stars 
during the main burst. 

In conclusion, for the most massive spheroids, the interplay 
between the peak in the merger rate and the subsequent AGN 
quenching of the star formation act together in such a way that 
most of the star formation process and the galactic assembly oc- 
cur at roughly the same time and the same place, thus mimicking 
a sort of "monolithic" behaviour. In other words, even though the 
duration of the star formation would lead to quite low [a/Fe] ra- 
tios from the point of view of a pure chemical evolution model, 
the fact that it had happened in several sub-units makes the final 
[a/Fe] ratios higher and in better agreement with observations. 
Intermediate and small objects, instead, do not have a quenching 
mechanism acting directly at their scales which can self -regulate 
the duration of the star formation. Therefore they end up having 
either too high or too low [a/Fe] ratios. 

In the remainder to this section we discuss the main ingre- 
dients that might be modified in order to improve the predicted 
[a/Fe]-mass relation. 

6.1. AGN feedback 

We start presenting the key factor that permitted substantial im- 
provement. If we switch off the SMBH heating of the intra-halo 
gas, we predict that the more massive galaxies are younger than 
the less massive ones - with a typical age of 6 Gyr (at variance 
with the observational results) - and that they are strongly a- 
depleted. As can be seen in fig.[T3](top panel), we basically con- 
firm with a more self-consistent model the results that Thomas 
(1999) obtained: the natural prediction of a bottom-up hierarchi- 
cal growth of the galaxies leads to a a/Fe -mass relation which 
has a negative slope. 

Therefore some kind of SF quenching is needed. The recipes 
commonly adopted in the literature (see Introduction) and in the 
present work assume that only the most massive galaxies are af- 
fected by SMBH feedback. This helps to predict higher [a/Fe] 
ratios and relaxes, but does not solve yet, the problem at the high 
mass. As discussed above, however, little attention is paid to 
what becomes of the less massive objects which are the building 
blocks of the relation and responsible for the bulk of the scatter. 
Suppression of star formation at low and intermediate masses 
is required in the models. In fact, recent observations show that 
AGN feedback appears to be present also in this mass range of 
the early-type galaxy population (Schawinski et al 2007) 

The bottom panel of fig. [T3] show us that the predicted MMR 
slope is in better agreement with the data, although offset down- 
wards in line with the above discussion (see Sec. 5) and PMOS's 
analysis. 

6.2. Stellar yields 

A change in the stellar yields will introduce a systematic offset 
of a few tenths of a dex in the model predictions (see Thomas et 
al., 1999, PM04), hence it might let all the massive galaxies lie 
within the observational boundaries for a suitable choice of the 
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Fig. 13. The o-ZFe-mass (top panel) and the MMR (bottom 
panel) relations as predicted by GallCS for the whole sample 
of ellipticals (black points). In this run the AGN feedback has 
been switched off. The thick solid lines encompass the Icr-region 
(hatched area) around the mean trend reported by Thomas et al. 
(2008). Note that in this case we have fewer ellipticals than in 
fig. [3] simply because gas is allowed to cool onto a disc at the 
centre of massive DM haloes, which leads us to classify more 
central galaxies as disc-dominated spirals. 



stellar nucleosynthesis is donq3- However, being only an offset, 
this change will not modify the slope of the predicted ajFe- 
mass relation and will exacerbate the problems for the low mass 
galaxies. But most importantly, the successful calibration of our 
model with element ratios observed in Milky Way stars does not 
allow significant modifications of the underlying stellar yields. 

6.3. IMF 

We do not test other IMFs, since a Salpeter IMF is a valid as- 
sumption for explaining most of the properties of early-type 
galaxies (see Renzini 2005). Furthermore, it has been akeady 
shown by Nagashima et al. (2005) that a change in the IMF is 
not enough (in the context of hierarchical mergers) in reproduc- 
ing the correct a/Fe-mass relation. We expect that a flattening 



^ For instance by: i) extending the upper mass limit of the IMF to 
100 Mo, ii) neglecting the ejecta of the stars in the mass range 8-11 Mq 
(whose O/Fe is slightly sub-solar); iii) by setting the SNIa rate to the 
lowest value permitted by observations. 



of the IMF regardless of the galaxy mass leads to a overall shift 
towards higher values of the predicted [a/Fe] ratio, but it does 
not affect the slope of the a/Fe-mass relation unless one finds 
a good reason to make the IMF flatter as the galactic mass in- 
crease. We will investigate this possibility in future work. This 
solution appears contrived, however, because it implies that the 
single building block should know in advance its destiny in order 
to self-assign a suitable IMF. 

6.4. Feedback from SNe and its efficiency 

In Eq. (7) e"' represents the efficiency of mass-loading during 
the triggering of a galactic wind by SNII explosions. Decreasing 
e produces more feedback, heating more cold gas, ejecting more 
hot gas from halos, and thus reducing the amount of gas that 
can potentially form stars. In this case we find that the predicted 
stellar masses are smaller than the fiducial case. The galaxies 
look slightly more a-enhanced as expected since the SF process 
is strongly-disfavored by the SNe explosions. However, this is 
not a viable solution for the a/Fe-mass relation problem, since 
the high mass-loading also implies a very low metal content in 
the stars. The predicted MMR is offset downwards by at least 0.5 
dex from the observational one. On the other hand, if we switch 
the SN feedback off, we tend to slightly worsen the a/Fe-mass 
relation, whereas the agreement for the MMR improves. 

Following the above line of thoughts, we further modify 
GallCS by introducing SNIa contribution into Eq. |7] Since a 
elements and Fe are still ejected at the same rate, this change 
has the same effects of decreasing 43- A 0. 1 dex increase in the 
final [a/Fe] ratios can be achieved when a differential wind is 
invoked, namely if we assume that twice more Fe than O can be 
ejected in the hot phase due to SNIa explosions. Again, given 
the nature of such a mechanism, neither the slope of the pre- 
dicted a/Fe-mass relation can be steepened nor its scatter re- 
duced. Further investigation will tell us if a change in the SNe 
feedback, namely by allowing them to quench the star formation 
as in monolithic models (e.g. Pipino et al., 2008b), might be the 
required galactic-scale source of feedback. 



6.5. Star formation efficiency 

We identify the SF recipe as one of the prescriptions where 
one can improve upon. In fact, Pipino et al. (2008b) started 
from the heuristic approach of PM04, who required the SF ef- 
ficiency to increase as a function of galactic mass, and showed 
that the a/Fe-mass relation can be explained by implementing 
a physically motivated value for the SF efficiency. To explain 
the higher star formation efficiency in the most massive galax- 
ies, they appeal to massive black holes-triggered SF: a short 
(10'' - 10^ yr) super- Eddington phase can provide the acceler- 
ated triggering of associated star formation. The SMBH grows 
mostly in the initial super-Eddington phase while most of the 
spheroid stars grow during the succeeding Eddington phase, un- 
til the SN-driven wind quenches SF. According to Pipino et al. 
(2008b) models, the galaxy is fully assembled on a time-scale 
of 0.3-0.5 Gyr. This time-scale is long enough, however, to al- 
low the SMBH to complete its growth in order to reproduce the 
MagoiTian relation. The fact that GallCS akeady turns gas into 
stars at the maximum possible rate during the merger-induced 
star-burst phase and its quite low mass/space resolution hamper 
us from a direct implementation of the above recipe. Moreover, 



^ With the obvious difference that we can have ejection of matter also 
when the SF is zero, because of the nature of SNIa progenitors. 
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the fact that stars born the disc and can be transferred to the 
bulge of the same galaxy because of instabiUties, is a possibil- 
ity not taken into account in Silk (2005). This scenario will be 
tested in the forthcoming version of GaUCS. 



7. Conclusions 

We implemented a detailed treatment for the chemical evolu- 
tion of H, He, O and Fe in GallCS, a semi-analytical model for 
galaxy formation which successfully reproduces basic low- and 
high-redshift galaxy properties. The contribution of supemovae 
(both type la and II) as well as low- and intermediate-mass stars 
to chemical feedback are taken into account. The model pre- 
dictions are compared to the most recent observational results 
by Thomas et al. (2008). We find that this chemically improved 
GallCS does not produce the observed mass- and cr-[a/Fe] rela- 
tions. The slope is too shallow and scatter too large, in particular 
in the low and intermediate mass range. The model shows sig- 
nificant improvement at the highest masses and velocity disper- 
sions, where the predicted [a/Fe] ratios are now marginally con- 
sistent with observed values. Moreover, an excess of low-mass 
ellipticals with too high a [a/Fe] ratio is predicted. We show that 
this result comes from the implementation of AGN (plus halo) 
quenching of the star formation in massive haloes. 

A thorough exploration of the parameter space shows that 
the failure of reproducing the mass- and o--[afFe] relations can 
partly be attributed to the way in which star formation and feed- 
back are currently modelled. The merger process is responsible 
for a part of the scatter. We suggest that the next generation of 
semi-analytical model should feature feedback (either stellar of 
from AGN) mechanisms linked to single galaxies and not only 
to the halo, especially in the low and intermediate mass range. 

Furthermore, a drawback of the is the fact that the MMR can- 
not be fit simultaneously. Both effects can be explained by the 
fact that the model is still lacking a sort of monolithic formation 
for all its spheroids, which is needed in order to reproduce the 
a/Fe-mass relation and the MMR at the same time. The scatter 
is somehow intrinsic to the merger history, thus calling for fur- 
ther modification of the baryons behaviour with respect to the 
CDM. In other words we envisage a lack of a self-regulating 
mechanisms which acts on a galactic scale and counterbalances 
to some extent the random nature of the merger trees. 

As expected from chemical evolution studies, is the shape of 
the SFH which sets the final [a/Fe]: a galaxy with a shorter du- 
ration of the SFH (summed over all the progenitors) will have a 
higher [a/Fe] than a galaxy with a longer one, even if the latter 
had less mergers. Moreover the [a/Fe] achieved by the galax- 
ies are in general 0.1-0.3 higher than what expected by feeding 
the integral SFH in a pure chemical evolution model. This hap- 
pens because in GallCS galaxies do not evolve as closed boxes. 
They instead exchange metals with the surrounding hot halo, un- 
dergo dry-mergers, as well as stars can be created in discs and 
the moved to bulges because of mergers or instabilities. 

In order to understand such a difference and to find viable so- 
lutions we tested the effect of several model parameters. Among 
those, we emphasise that an increase in the star formation ef- 
ficiency and Fe-enhanced winds driven by the SNla activity 
might play a role in removing galaxies with too a low [a/Fe] 
ratio. However, given the way they act on the galaxy evolution, 
the cannot be effective in either making the slope of the pre- 
dicted or/Fe-mass relation steeper or in reducing its scatter. In 
particular, it seems hard to remove the low-mass galaxy too a- 
enhanced. 
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